TeV 7-rays and neutrinos from photo-disintegration of nuclei in Cygnus OB2 
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TeV 7-rays may provide significant information about liigh energy astrophysical accelerators. Such 
7-rays can result from the photo-de-excitation of PeV nuclei after their parents have undergone 
photo-disintegration in an environment of ultraviolet photons. This process is proposed as a can- 
didate explanation of the recently discovered HEGRA source at the edge of the Cygnus OB2 as- 
sociation. The Lyman-a background is provided by the rich O and B stellar environment. It is 
found that (1) the HEGRA flux can be obtained if there is efficient acceleration at the source of 
lower energy nuclei; (2) the requirement that the Lorentz-boosted ultraviolet photons can excite 
the Giant Dipole Resonance implies a strong suppression of the 7-ray spectrum compared to an 
E~'^ behavior at energies < 1 TeV (some of these energies will be probed by the upcoming GLAST 
mission); (3) a TeV neutrino counterpart from neutron decay following helium photo-disintegration 
will be observed at IceCube only if a major proportion of the kinetic energy budget of the Cygnus 
OB2 association is expended in accelerating nuclei. 



I. INTRODUCTION 

There are two well-known mechanisms for generating 
TeV 7-rays in astrophysical sources [Ij] . The first involves 
purely electromagnetic (EM) processes, including syn- 
chrotron emission and inverse Compton scattering. The 
second may be termed hadronic, in which the 7-rays orig- 
inate in tt" decays. The latter can in turn be traced to 
either pp or p7 collisions. 

The EM processes originate in the acceleration of elec- 
trons; as a result, the TeV 7-rays must be accompanied 
by an X-ray counterpart. In addition, the conflicting re- 
quirements on the magnetic field (large field for acceler- 
ation, small field for limiting synchrotron cooling) limits 
the energy of the emitted photons E.^ < 10 TeV. 

The hadronic pj mode is characterized by a large 
threshold for vr*^ production, and thus it is favored only in 
very hot photon environments, or in the presence of very 
energetic proton beams. In addition, the presence of the 
high threshold allows acceleration free of scattering losses 
as long as the Hillas criterion holds. In the pp hadronic 
mode, threshold effects are insignificant. Because of this, 
it is generally assumed that tt" production occurs in a 
region distinct from the site of the primary acceleration. 
However, one can set conditions relating the interaction 
length of the nucleons, the neutron decay lifetime and the 
source confinement radius to permit simultaneous emis- 
sion of neutrons, 7-rays and neutrinos In the pj mode 
the spectrum of 7-rays follows that of the parent proton 
population in the energy region above threshold. In the 
pp mode, the resulting 7-ray spectrum is broad, refiect- 
ing the presence of a quasi- Feynman plateau which spans 
the entire rapidity space. The photon spectrum will de- 
viate from the proton spectrum if the Feynman plateau 
is not fiat in rapidity space (^], as suggested by Teva- 
tron data Finally, for both the pj and pp modes, the 
decay of the charged pions yield a neutrino counterpart 



with energy and intensity similar to the photons. 

In this paper, we discus in detail a third dynamic which 
can lead to TeV 7-rays: the photo-disintegration of nu- 
clei at the source, followed by the photo-de-excitation of 
the daughter nuclei In order to generate TeV 7-rays 
as a result of emission of MeV 7-rays in the rest frame of 
the de-exciting nucleus, the Lorentz factor of the boosted 
nucleus must be ~ 10®. For this boost factor, excitation 
via the Giant Dipole Resonance (~ 10 MeV — 30 MeV in 
the nucleus rest frame) is obtained with ambient photons 
with energies in the far ultraviolet (usually defined as 1- 
20 eV). Photons of these energies are expected from the 
Lyman a emissions from hot stars. This process clearly 
reduces the threshold energy requirement relative to p^. 
The important role played by the Giant Dipole Reso- 
nance (GDR) in the photo-disintegration effectively sup- 
presses the 7-ray spectrum below 1 TeV. 

In recent decades, the Cygnus Spiral Arm has been a 
site of 7-ray signal candidates which tend to come and 
go. In the energy band 2 x 10® < E^/GeV < 2 x 10^, 
data collected by the Kiel air shower experiment Q show 
a 4.4cr excess of events in the direction of the binary 
system Cygnus X-3, with the typical 4.8 hr modulation 
previously observed in the MeV [1| and TeV regions. 
This result, including the X-ray binary period, was later 
confirmed by the Haverah Park data |10j], with the ad- 
ditional observation of an abrupt steepening of the spec- 
trum a,t > 2x 10^ GeV. The distance (« 10 kpc) to 
Cygnus X-3 slightly exceeds the minimum path length for 
cosmic microwave background (CMB) absorption in the 
energy band of the Kiel experiment [lli|. However, since 
absorptive effects decrease for energies beyond 10*" GeV 
(see Fig. [T]), the observed steepening must be traced to a 
softening of the injection spectrum. 

At a much higher energy > 5 x 10* GeV, the analysis 
of the cumulative Fly's Eye data also revealed an excess 
of events from the direction of Cygnus X-3, with chance 
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probability of 6.5 x 10~^ A 3.5a excess from this 

direction and in the same energy region has also been 
observed at the Akeno air shower array [13] ■ The inferred 
signal fluxes are consistent at the Itr level. However, the 
evidence for the 4.8 hr modulation is termed as "weak" by 
the Fly's Eye Collaboration, and is absent in the Akeno 
data. Thus, one may infer that the high energy signals 
originate in a different source within the angular field of 
view of Cygnus X-3. 

More recent data from the CASA-MIA [ll| and the 
HEGRA [T3| experiments place restrictive bounds on 
steady state fluxes from the Cygnus region. These can 
be seen in Fig. [TJ and suggest that the earlier reported 
fluxes do not reflect current steady state activity. At this 
point, it is critical to note that only the HEGRA exper- 
iment has the angular resolution to place an upper limit 
on the steady state flux from Cygnus X-3, and at the 
same time observe significant activity (at the 7a level) 
in the TeV region from an unidentified source (with no 
optical or X-ray counterparts) which is within 0.5° of 
the X-ray binary An excess at the 3.3cr level from 
the direction of this unidentified source is also present in 
the Whipple data [1^]. The strength and specificity of 
the source, and the distinct absence of an X-ray coun- 
terpart, make this source a good candidate for probing 
the nucleus photo-disintegration/de-excitation model for 
producing TeV 7-rays. 



Especially intriguing is the possible association of the 
TeV HEGRA source with Cygnus 0B2, a cluster of sev- 
eral thousands of young hot OB stars. At a relatively 
small distance (« 5000 light years) to Earth, this is 
the largest massive Galactic stellar association. Cosmic 
ray nuclei are expected to be trapped and accelerated 
through turbulent motions and collective effects of star 
winds. In this paper we present a model for explaining 
the HEGRA observations, in which the trapped high en- 
ergy nuclei undergo stripping on the starlight background 
and their surviving fragments emit 7-rays in transition 
to their ground states. A short version highlighting the 
salient features of the model has been issued as a compan- 
ion paper The outline here is as follows: in Sec. |TT] 
we provide a detailed description of the salient charac- 
teristics of the Cygnus 0B2 association, including stel- 
lar counts, lifetime, size and energy considerations. Sec- 
tion lHII contains a description of the HEGRA TeV source, 
and summarize the different proposed explanations. The 
model is presented in Sec. llVi including a calculation of 
the TeV photon flux, following a discussion of the accom- 
panying neutrino flux. In Sec. |V] we compare the TeV 
7-ray and neutrino yields in the model with those corre- 
sponding to the hadronic modes. In Sec. I VII we examine 
the possibility whether our model can accommodate the 
emission of ultra-high energy cosmic rays. Conclusions 
are collected in Sec. I VIII 
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FIG. 1: Integrated flux from the binary system Cygnus 
X-3 as reported by the different collaborations. The dot- 
ted line is a best single power law flt, Fj{> Ej) = 3.3 x 
10"** {Ej/GeV)~°-^^ cm-2 s~\ to the Cygnus X-3 data: 
A ^, ■ 'c?], • 1?], and T 0. The accompaning solid and 
dash-dotted lines take into account attenuation on the CMB, 
with the latter also including a cutoff at the source. Also 
shown are the integrated fluxes of neutral particles from the 
direction of Cyg nus X-3 reported by the Akeno (□ and 
Fly's Eye (★ [12]) collaborations. The cross-hatched bands in- 
dicate the 90% CL upper limit on steady state flux from the 
direction of Cygnus X-3 as observed by the CASA-MIA ^ 
and HEGRA experiments. 



II. CYGNUS OB2 IN A NUTSHELL 

The Cygnus 0B2 (VI Cygni) association is one of the 
most massive associations, with some of the most lumi- 
nous stars, in our Galaxy. It was first noticed by Miinch 
and Morgan TSj who, during the course of a classification 
of blue giants (OB stars), found eleven of these objects 
in that region. Since the pioneering spectroscopy and 
photometry of Johnson and Morgan [19] and subsequent 
works [20| . this region has been known to harbor a large 
population of massive and early type stars, which have 
been found to be highly reddened. The first comprehen- 
sive study of this stellar association 21, 22, 23] identified 
a few hundred OB stars as possible members based on 
photographic photometry. Other photometric and spec- 
troscopic studies [2^ were carried out during following 
years confirming and extending the first results. 

Cygnus 0B2 is located at galactic coordinates (/, 6) ~ 
(80°, 1°), behind the Great Cygnus Rift. The study per- 
formed in Ref. (2^ . inferred an elliptical shape with ma- 
jor and minor axes of 48' and 28', respectively, and esti- 
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mated more than 3000 stars with at least 300 of OB spec- 
tral type, which resulted in a total mass for the associa- 
tion of (0.6 — 2.7) X 1O^M0. However, due to the extreme 
redenning in and around this region, which hampers the 
detection of even OB stars, the observed morphology was 
rather an artifact created by the particular extinction 
pattern of the field. A more recent work based on a 
statistical study of point sources revealed by the Two Mi- 
cron All Sky Survey (2MASS) in the near-infrared, shows 
that this association is much larger and richer than previ- 
ously thought. The resulting stellar distribution reveals 
a rather regular and almost circular density profile with 
the center located at ia,S) = (20''33"1G", -f41°12') and 
with a pronounced maximum slightly offset at {a, 6) — 
(20^33'"10^-h41°15.7'). Stars counts show that 50% of 
the members are located within a radius of 21', and 90% 
within a radius of 45' around the center, merging with 
the field stars at a radius of ~ 1°. The central stel- 
lar density reaches 4.5 stars/arcmin^ above the field star 
density, and drops to 50% at a radius of 13'. By inte- 
grating the radial density profile, after substraction of 
the field star density, the total number of OB stars was 
found to be 2600 ± 400 [1^ , much larger than the earlier 
estimate of ~300 Furthermore, the total number 

of 0-type stars was inferred to be 120 ± 20 Nev- 
ertheless, this number must be cast as an upper limit, 
as the 2MASS data do not allow a precise spectral de- 
termination and hence do not account for possible evolu- 
tionary effects, i.e., the quoted number reflects the initial 
number of very massive (> 20Mq) stars, although some 
may have already evolved. This also agrees with more 
recent estimates [1^ [l^l and suggests a total mass for 
the association of (4 — 10) x 10^ M0, where the primary 
uncertainty comes from the unknown lower mass cut-off. 
Using the radial density profile, a central mass density of 
(40 - 15O)M0/pc3 is determined [25|. 



Early distance determinations revealed the proximity 
of this region, varying from 1.5 kpc fl§\ to 2.1 kpc [2^ . 
More recent estimates have set this distance to d ~ 
1.7 kpc [H, Hi]. At such distance, the inner 21', with 
half the total number of objects, results in a physical ra- 
dius of 10 pc, with i?out ^ 30 pc being the radius 
of the association. Yet, a revision of the effective temper- 
ature scale of 0-type stars yields a closer distance of '--^ 1.5 
kpc [l^l , when the age of the association is taken into ac- 
count in fitting the late-0 and early-B dwarfs to model 
isochrones. It is from the superposition of the isochrones 
(calculated using the theoretical evolutionary tracks of 
the model of Ref. [30*1) on the Hertzsprung- Russell dia- 
grams, that the age of the association has been estimated 
to lie between 1-4 Myr [3l[, compatible with other esti- 
mates [13, [2^ [13, [131 . This range reflects the dispersion 
of the upper main sequence and agrees with the fact that 
the observed large number of 0-type stars implies that 
the association should be younger than '^5 Myr, because 
in the case of coeval star formation, the number of this 
type of stars decreases rapidly [ll|. In addition, the fact 
that there are some Wolf-Rayet stars within Cygnus 0B2 



suggests an age larger than ~2 Myr. On the other hand, 
the non-detection of any supernova remnant in Cygnus 
0B2 [1^ points to an association younger than ~4 Myr. 
Nevertheless, the presence of a few slightly evolved su- 
pergiants (an 03 If* star, three Wolf-Rayet stars and 
two LBV candidates) suggests that the star formation in 
the asociation was not strictly coeval and there arc even 
indications of ongoing star formation [28,, J9., jiS,, j35, 36J. 
This fact would narrow down the quoted age boundaries. 

Early estimates of the wind mechanical luminosity for 
the Cygnus 0B2 association gave L^j — 10'^* erg s~^ [13]. 
However, the improved data on the population of the re- 
gion required a revision of this estimate. By using a de- 
tailed analysis of radiatively wind models for hot stars, 
it has been shown than at an early stage (< 2 Myr), 
the wind mechanical luminosity is maintained approxi- 
mately constant [13], being ^ 2 x 10"^"' erg s~^ per solar 
mass, which, for the association, yields L^, ~ (1 — 2) x 
j^q39 ej-gs~^ i3S]. After the flrst ~ 2 Myr, the luminosity 
increases further as Wolf-Rayet stars appear. This esti- 
mates assumes stars of solar metallicity and a Salpeter 
initial mass function (IMF) slope. There have been dif- 
ferent calculations for the IMF slopes [H, [H, [U [li] 
which however do not agree with each other, although 
lie in the interval F ~ -(1.0 - 1.6) (F5 = -1.35 being 
the Salpeter slope). Nevertheless, these discrepancies are 
due to systematic uncertainties and hence, a reasonable 
approximation is to assume a uniform canonical Salpeter 
slope. Should the metallicity be close to the solar value, 
this implies ,37.] that the estimate of the wind mechani- 
cal luminosity is probably correct within a factor of ^ 2. 
Other values for the metallicity would modify this pre- 
diction by at most a factor < 10. On the other hand, the 
total Lyman continuum luminosity of the stars has been 
estimated to be ~ lO^^photons s~^ [Hf. 

Recently, there have been observations of TeV 7-rays 
from the northeast boundary of this association (see 
Fig. [2]) [15I] . These observations will be discussed in the 
following section, and will constitute the focus of this 
paper. 



III. TeV J2032+4130 

The HEGRA system of Imaging Atmospheric 
Cerenkov Telescopes (HEGRA lACT-system) consisted 
of five identical Cerenkov telescopes (each with 8.5 m^ 
mirror area) and employed a stereoscopic technique 
achieving an angular and energy resolution better 
than 0.1° and 15%, respectively, for 7-rays on an 
event-by-event basis with energies from 0.5 TeV to ^ 
50 TeV [H]. 

The observation of the Cygnus region by the HEGRA 
lACT-system has allowed the serendipitous discovery of 
a TeV source in the outskirts of the core of Cygnus 
0B2 [13]. The analysis of the total 278.3 hours of ob- 
servations performed in two periods from 1999 to 2002 
(120.5 hours from 1999 to 2001 ^ and 157.8 hours dur- 
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ing 2002 [l5[) has revealed the presence of a steady (and 
possibly extended) TeV source, with hard injection spec- 
trum. Interestingly, there have been earlier claims of a 
multi-TeV excess in this region [41*1. 

The excess significance of the TeV source is 7.1a and it 
appears extended at more than Aa level with a morphol- 
ogy which is suitably described by a Gaussian profile. 
The source is termed TeV J2032-I-4130 after the position 
of the center of gravity. Its extension (Gaussian Icr ra- 



dius) is 6.2'(±1.2^j^j ± O.g^yJ, which at a distance of 1.7 
kpc results in r ~ 3.07(±0.59stat ± 0.45sys) pc. 

Three different types of intrinsic morphology were 
tested [15] (disc, volume and surface), but the data can- 
not discriminate between them. The energy spectrum 
determination yielded a pure power-law fit with a hard 
photon index, showing no indication for an exponential 
cut-off, given by [lH| 
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which implies a flux above 1 TeV given by 

F^{E^ > ITeV) = 6.9 (±1.8stat) x lO-^^cin-^s-\ (2) 

These results imply a luminosity of ~ 10^^ erg/s above 
1 TeV, which is well within the kinetic energy budget of 
Cygnus 0B2, and indeed, also within that of a number 
of notable member stars [29l . [43| . 

As can be seen in Fig. [21 this signal is located at 
the edge of the error circle of the EGRET source 3EG 
J2033-I-4118 and within the ^10 pc-radius core circle 
of the Cygnus 0B2 association. The EGRET observa- 
tions [45!| lie in the energy range below 10 GeV, and 
closely display a spectrum cx E~'^. It will be of important 
significance for our discussion that when extrapolated to 
the TeV range according to the E~'^ behavior, the fiux 
exceeds the HEGRA fiux by several orders of magnitude. 
Hence we surmise that these observations constitute two 
different sources. 

An additional set of observations performed dur- 
ing 1989-90 by the Whipple Observatory atmospheric 
Cerenkov imaging telescope 16] has been recently rean- 
alyzed in the light of the HEGRA data. These confirm 
an excess in the same direction as J2032-f4130, although 
with considerably larger fiux, above a peak energy en- 
ergy response of 0.6 TeV. The statistical significance of 
the signal is only 10% smaller with selection of events 
above 1.2 TeV. However, the large differences between 
the fiux levels cannot be explained as errors in estima- 
tion of the sensitivity of the experiments since they have 
been calibrated by the simultaneous observations of other 
TeV sources. 

So far no clear counterparts at other wavelengths have 
been identified, and moreover, the observed spectrum is 
not easily accommodated with synchrotron radiation by 
electrons [13] • The difficulty to accommodate the spec- 
trum by conventional electromagnetic mechanisms has 
been exacerbated by the failure of CHANDRA and VLA 
to detect X-rays or radiowaves signaling acceleration of 
any electrons [441 ]. 

Nevertheless, a leptonic origin is not yet excluded, es- 



pecially if the TeV source is actually located in the vicin- 
ity of Cygnus X-3. The suggested model [40| proposes a 
jet-driven termination shock at the boundary where the 
relativistic jet meets the interstellar medium. This re- 
sults in synchroton and TeV inverse-Compton emission 
from the accelerated electrons [i^ , giving rise to the ob- 
served signal. Such a jet could emanate either from an 
as-yet undiscovered microquasar or from Cygnus X-3. In- 
terestingly, the TeV signal aligns well with the northern 
error cone of the bi-lobal jet of Cygnus X-3 [13] ■ For the 
observed angular separations, this would place the TeV 
source at a distance of 10 kpc from Earth, and around 
70 pc from Cygnus X-3. 

The fact that there is no catalogued X-ray sources 
within the 2a error circle of the TeV source would dis- 
favor a leptonic origin at Cygnus 0B2. We now briefly 
discuss some models in which accelerated protons and nu- 
clei interacting with a local dense gas cloud [3] or with 
stellar winds [48| , produce photons from 7r° decay. In the 
case of gas cloud interactions [ij , an interstellar density 
of ntot 30 cm^'^ and an injection efficiency (the ratio 
of cosmic ray energy to the kinetic energy of the acceler- 
ating winds) of 0.08% reproduces well the observed TeV 
spectrum. However, such a large density has been called 
into question [48] because a re-assessment [i^ of the CO- 
H2 conversion factor for this type of environment (used 
for calibration) would imply a considerably lower den- 
sity, ntot ~ 0.1 cm~'^. With such a density, a much larger 
injection efficiency 25% would be required. A differ- 
ent model in terms of cosmic ray illumination of stellar 
winds has been suggested to explain the signal through 
hadronic interactions in the innermost parts of the winds 
of massive OB stars [i^. However, a more recent de- 
tailed study of the stellar winds [13], including pressure 
effects and a more accurate estimate of the (lower) mass- 
loss rate for OB stars, resulted in a significant reduction 
in the photon emission rate. Thus, if only the currently 
known stars exist, the predicted flux from this model is 
too low to comfortably explain the TeV source. 

Another suggested possibility has been that of a point- 
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Cygnus 0B2 Field: HEGRA CT-System 
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FIG. 2: Skymap of correlated event excess significance (a) 
from all HEGRA lACT-System data (3° x 3° FoV) centered 
on the TeV source J2032+4130. Nearby objects are also 
shown: 95% contours for 3 EGRET sources (indicated by 
the ovals), their possible X-ray associated counterparts (as 
given in Ref. 53 )> ^-^d Cygnus X-3. The center of gravity 
with statistical errors and intrinsic size (standard deviation 
of a 2-dim Gaussian, (Jsrc) are indicated by the white cross 
and white circle, respectively. The TeV source, J2032+4130, 
is positioned at the edge of the error circle of the EGRET 
source 3EG J2033+4118, and within the core circle of the 
extremely dense OB stellar association Cygnus OB2 [isl ]. 



like origin of the TeV signal related to an unusual tran- 
sient X-ray source which lies 7' from the center of gravity 
of TeV J2032-h4130 jSll]- However, the position and vari- 
ability of this X-ray source and mainly, the confirmed 
extension of the TeV source makes this a remote alterna- 
tive, although it remains possible that several point-like 
TeV sources could be masquerading as a single extended 
source. 

The TeV signal could also be interpreted as originated 
in the wind nebulae of an as- yet undetected recent pulsar 
(similar to Vela-type pulsars) , as a result of hadronic or 
leptonic interactions |52l |. This pulsar must have been 
formed ~ 10'* years ago in a core collapse of one of the 
massive stars in the Cygnus 0B2 association. However, 
the non-observation of any supernova remnant in the re- 
gion, as well as the greatly lowered estimates of the ambi- 
ent gas density, seem to pose a serious challenge for this 
model. 

In summary, there is no single compelling explanation 
for the various characteristics of TeV J2032-I-4130 [s^. 
In the next section, we present a new mechanism for ex- 
plaining the HEGRA data. 



IV. THE MODEL 

There are two processes by which the nucleus may 
produce 7-rays, photo-disintegration of the nucleus by 
ambient photons followed by de-excitation of a daughter 
nucleus, and photo-pion production followed by decay 
of the neutral pions. The two processes have different 
thresholds and different signatures. In this section, we 
are concerned with the former mechanism, emphasized 
more than a decade ago by Moskalenko and collabora- 
tors 
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but largely ignored by the rest of the 7-ray 
community. In the following section we will compare and 
contrast the two processes. 



A. Nucleus Photo-Disintegration 

The interaction between photons and high energy nu- 
clei results in the emission of nucleons. The relevant 
photonuclear interaction process for the relevant ener- 
gies have been studied from the point of view of the 
collective and the shell models [5a]. It is shown that 
collective nuclear states dominate the interaction, with 
low angular momentum modes preferred. In the energy 
region which extends from threshold for single-nucleon 
emission - 10 MeV up to 30 MeV the GDR dom- 
inates. The GDR typically de-excites by the statisti- 
cal emission of a single nucleon. Above the GDR re- 
gion, and up to the photo-pion production threshold at 
E^Yi — "Ztt (1 + m7r/2mAr) ~ 145 MeV, the non-resonant 
processes provide a much smaller cross section with a 
relatively flat dependence on energy. 

The photo-disintegration rate for a highly relativistic 
nucleus with energy E = jAmN (where 7 is the Lorentz 
factor) propagating through an isotropic photon back- 
ground with energy e and spectrum n(e), normalized so 
that the total number of photons in a box is / n(e)de, 
is M 



c 1 



2 2 
7^6^ 



de 



/ e'crA{e')de, 

JQ 



(3) 



where o'^(e') is the cross section for photo-disintegration 
of a nucleus of mass A by a photon of energy e' in the 
rest frame of the nucleus. 

The cross section for all the different nuclear species 
has been obtained through a direct fit to data [STj . 
For medium and heavy nuclei {A > 30) the total pho- 
ton absorption cross section can be approximated by 
a dipole form (sometimes called a "Breit-Wigner" or 
"Lorentzian" ) 



cr^(e') = ctq 



(ef,2-e'2)2 + £'2r2 ' 



(4) 



where T is the width, ep is the central value of the GDR 
energy band, and ctq is the normalization. We have found 
that for the considerations in the present work, the cross 
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section can be safely approximated by the single pole of 
the Narrow- Width Approximation, 



r 

crA(e') = TT do - 6{e' - Eq) , 



(5) 



where ao/A = 1.45 x 10"27(,j^2^ T = 8 MeV, and e'n = 
42.65yl-°-2i (0.925^2.433) ^^y^ for A > 4 (A < 4) [H, 
[sgj . Inserting Eq. Q into Eq. ^ we obtain 



„ TT (To fn r f°° de , X ^ ,„ , X 
47 Jo e 



472 

TT gp gp r 
472 



0/27 



rfe 



(6) 



As a test of our approximation, we first compute the 
disintegration rate for a nucleus passing through a region 
where 



n(£)=n?E(6) = (e/7r)^ 



c/T 



1 



(7) 



corresponding to a Bose-Einstein distribution with tem- 
perature T. The result is 



R 



BE 



In f 1 - e-^o/27T 



(8) 



We then verify that for '^^Fe and for the CMB tempera- 
ture (T = 2.3 X 10~^ eV), this solution agrees to within 
20% with the parametrization given in [60|. The latter 
was derived using the full form of the cross section pre- 
sented in [H^l • (A similar result in the context of photo- 
disintegration of nuclei at the Galactic center has been 
obtained in [sij.) 

In the companion paper [l7| it was argued that the rate 
can be written as a function w'^ \ ln(l — e"*")! of a scaling 
variable w = eQ/2jT, times a prefactor {ao T / e'q) / n . 
The peak rate therefore scales as the prefactor. Using 
the fits mentioned below Eq. (5), the prefactor in turn 
scales in A as A^-^^, with a small correction for A < 4 
nuclei such as helium. 



B. The Photon Population 

The ingredients necessary for calculating the total 
photo-disintegration rate from a given region of the OB 
association are (i) an ambient photon distribution in or- 
der to obtain the rate R\ on starlight per nucleon, and 
(ii) an initial population density of nuclei ua in the re- 
gion. In a qualitative manner, the iron (silicon) photo- 
disintegration rate in a region of radius Rm « 10 pc 



dN 
~dt 



Fc(Si) ^56(28) 



-(4/3)7ri?f„ 7iFc(si) 



-^56(28) 



(9) 



where npo(si) is the iron (silicon) population density in 
some energy bin. This population will be assumed to re- 
sult from continuous trapping of the diffuse cosmic ray 



flux by diffusion in a milligauss magnetic field 'g?] , acting 
preferably on heavy nuclei. Energies of O (PeV/nucleon) 
are achieved through re-acceleration in strong winds of 
the OB stars. We assume, and justify a posteriori^ 
that the nucleus population n.p(,(gi) does not significantly 
change during the time considered. It will be seen below 
that at most one nucleon will he stripped in the region of 
interest during the diffusion time within the association; 
in the calculation of the disintegration rate, this justifies 
inclusion of only the lowest order in the photo-nuclear 
interaction. 

The photon background will be assumed to result from 
the thermal emission of the stars in the region R. The av- 
erage density in the region R will reflect both the temper- 
atures To E^nd Tb due to emission from O and B stars, re- 
spectively, and the dilution resulting from inverse square 
law considerations. Specifically, for a region with A^o O 
stars and A^b B stars, the photon density is 



'n^^ie)NoRl 



i22 



Nb Rl 



(10) 



where -Ro(B) is the radius of the O (B) stars, the factor 
9/4 emerges when averaging the inverse square distance 
of an observer from uniformly distributed sources in a 
region R . We take the O and B star populations to 
be 5% and 95%, respectively, of the total OB population 
of 2600 in the association. At this point, we do not dif- 
ferentiate between the 3 pc cell HEGRA hot spot and 
the rest of the cells in the association. In order to take 
into account the larger density of stars in the region of 
interest, we consider all the stars to lie within the central 
region, R = R^ ~ 10 pc. 

In Fig. [3] we show the dependence on the Lorentz factor 
of i?5g, i?2g, and R\ for the stellar ambiance described 
above (solid lines). The peaks of the three rate curves, 
while occurring at slightly different values of 7, are seen 
to scale as A^ '^ . This is in accord with our discussion of 
scaling below Eq. ([5]). Also in accord with our scaling 
formula, the positions in 7 of the peak rates should scale 
in A as 7poak ~ ep/^pcak ~ ^^°'2i. The three curves 
bear this out. We have taken for the O (supergiants [64] ) 
stars a surface temperature Tq = 40000 K, and radius 
Ro = 19 Rq; for the cooler B stars we assign Tb = 
18000 K and radius Rb = Rq For all three 

nuclei, the disintegration time exceeds the diffusion time 
(~ 10''^ yr [i^]) of the nucleus in the association. Thus, 
the a priori assumption of a lowest order calculation has 
been justified ^66*1. 

At first sight this model might seem to be a very crude 
approximation of the actual distribution of stars. We 
have also calculated the photo-disintegration rate with a 
more detailed model of the star distribution. We consider 
the thermal emission to come from the stars in the whole 
association, i?out 30 pc, and that half of them are 
uniformly distributed in the inner region, i?i„ 10 pc, 
and the other half are uniformly distributed in the outer 
shell, i. e. the density of stars in the inner region is 
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model we are proposing is stable. It is clear, however, 
that within the 3 pc HEGRA hot spot the concentration 
of stars would be above average, and thus hereafter we 
take as a fiducial value for R\ the one resulting from 



Eq. pO]) (more on this below). 

After the high energy nuclei interact with the photon 
field entering the GDR energy region, the nucleus is left 
in an excited state which will go o ver into an excited 
daughter state emitting 7-rays [5J]. Some early semi- 
quantitative statistical-model calculations for the pro- 
duction of 7-rays through the decay of the GDR in the 
^^Fe nucleus showed that the mean energy of the 7- 
spectrum is E'^^q ~ 2 — 4 MeV and the average multiplic- 
ity is 77.56 ~ 1 — 3 (68l |. Previous measurements showed 
that for the case of ^^O, the correspon ding values are 
E'^^Q - 5 - 7 MeV and riTi" ~ 0.3 - 0.5 Hence, in 

the observer system, these relativistic nuclei are a source 
of directional 7-rays with energy of the order ^ 7 MeV. 



FIG. 3; Photo-disintegration rate of ^''Fe, ^**Si, and ''He on 
the Cygnus 0B2 starlight. Solid (dashed) lines represent the 
simplified (more elaborate) model as described in the text. 



(_Rout/-Rin)"^ — 1 ^ 26 times that in the outer shell. The 
photo-disintegration is calculated to occur in a region of 
radius r ^ 3 pc in the outskirts of the inner part of the 
association, i?in, to model the size and position of the 
source of the HEGRA signal. This results in a photon 
density 



47 

^ ' 4 



^out 



(11) 



The factor 47/4 is a consequence of averaging the inverse 
square distance within this distribution for the density 
and the region where the reaction takes place [g^- In 
Fig. [3] we also show these results (dashed lines), which 
agree very well with the previous less elaborate model 
for the distribution of stars. This shows that even if we 
do not have a quantitative explanation for the anisotropy 
at the border, different global descriptions of the associa- 
tion point to the same photo-disintegration rate, and the 



C. TeV Gamma Ray Emission 

The low energy cutoff on R*^ seen in Fig. [3] will be 
mirrored in the resulting photon distribution. The E~'^ 
energy behavior of the various nuclear fluxes will not 
substantially affect this low energy feature. It is a ro- 
bust consequence of the model, discussed in and be- 
low Eq. p8)) . The energy behavior for photons in the 
1 — 10 TeV region of the HEGRA data is a complex con- 
volution of the energy distributions of the various nuclei 
participating in the photo-disintegration, with the rate 
factors appropriate to the eV photon density for the var- 
ious stellar populations. 

Let us define dRA/dE'^ as 



dRA 
dE' 



n(6) 



dt 



276 



d<T. 



7A 



dE'^ 



{e',E')de' , (12) 



where da^Ai^' , E'^)/dE'^ is the inclusive differential cross 
section for production of 7-rays from disintegration 
and E'^ is the energy of the emitted photon(s) in 
the rest frame of the nucleus. Assuming the same 
cosmic ray spectrum as above, the emissivity (num- 
ber/ volume/steradian) of 7-rays coming from nuclei de- 
excitation can be written as 



J 



Qf(i?,) = 5] j ^{E^)dE^ 



^ 5[E^ - 7S;(1 + cose^)] dE'^ 



(13) 



where E^ is the energy of the emitted 7-ray in the lab 
and is the 7-ray angle with respect to the direction of 



the excited nucleus. We write the nuclear flux as 



(14) 



8 



with Na a normalization constant, and Eq set to 1 TeV. 
Performing the angular integral with the delta-function 



constraint leads to 



J 



E 



rriN 
2 



N^i dE]\r 

"25 



duA , „ s dEN 
(En) 



E 



N 



^ dRA dE' 



dE^ 



(15) 



where Q is the Q-value of the de-excitation process. 
If we further approximate the 7-ray spectrum as be- 
ing monochromatic, with energy equal to its aver- 
age value (E'^j^), we can write da^A^ , E'^) / dE'^ — 

riAO'Aif') S{E!y — E'^^), where tTa is the mean 7-ray mul- 
tiplicity for a nucleus with atomic nurnber A. Hence, the 
emissivity can be approximated by [54| 



where 



£■^0 — 



'^E'^^Eq 
rriN 



2000S' 



fA 



(20) 



nArriN 



dEN 



En 



(16) 

The 7-ray emissivity is related to the differential flux at 
the observer's site (assuming there is no absorption) as 



dF-y 
dK 



47rd2 



(17) 



where Vdis is the volume of the source (disintegration) 
region and d is the distance to the observer. 

It is clear from Eq. (|16|) that if Ra is weakly depen- 
dent on En then the observed 7-ray flux will display 
the same power law behavior as the nuclei population. 
The HEGRA data show an approximate E^^ behavior 
for 1 TeV E^ < 1{) TeV, corresponding to a boost 
factor 10^ ^ 7 ^ 10^. Remarkably, as can be seen 
in Fig. [21 R\ varies by a factor of only 2 precisely in 
this region, far less than the 2 orders of magnitude of 
the primary nucleus flux. In order to evaluate Eq. (jl6p 
we approximate the behavior of R\ as roughly constant 
for i^min < E]S! < i^max, and zcro otherwise; where 
E^i^ ~ 10*^ GeV and E^e,^ - lO"^ GeV. Incorporat- 
ing Eq. ([H]) into Eq. p^ . we find that for a = 2 and 

^ 2E' . E„, 

E < 

• rtiM 



2 TeV (S;^/l MeV) 



A 



^E'^a 



R\ Na 



E„ 



Eo 



(18) 



independent of E^, whereas for Ej > 
2 TeV (E^/l MeV) 



2E' , E„ 



A 



^E'^A 



Na 



(19) 



The predicted constancy of the flux below ~ 2 TeV im- 
plies a strong suppression in this region relative to a 
flux extrapolated from HEGRA data to maintain an E^'^ 
behavior down to lower energies. For example, the up- 
coming GLAST mission will probe the gamma ray spec- 
trum from the Cygnus region in the range 20 MeV — 
300 GeV [7^. The flux predicted at ~ 100 GeV from an 
E~'^ extrapolation of the HEGRA data would render the 
source spectacularly visible in the GLAST observation, 
whereas the model here would predict a suppression by a 
factor of ~ (2/0.1)2 = 400 relative to this extrapolated 
flux. At 500 GeV the suppression would be still statisti- 
cally significant: E^'^x flux falls by a factor of ^ 16. 

The previous discussion has been qualitative. In Fig. [4] 
these features are display following a direct integration 
of Eq. (Uni), with choice of parameters to provide eyeball 
agreement with the data. One can see that the behavior 
of the fully integrated spectrum agrees in essential char- 
acteristics with the results of the preceding qualitative 
discussion. 

All in all, in the present model the photon flux follows 
the parent population of nuclei in the PeV/nucleon en- 
ergy region. Consistency then requires the nucleus spec- 
trum not to be significantly steeper than E~^. How does 
this comply with experimental data? 

As mentioned in Sec.[TTl the absence of supernova rem- 
nants in the association suggests that the nucleus popula- 
tion originates in the diffuse cosmic ray flux rather than 
in local supernova explosions. In assessing the consis- 
tency of the model we first obtain the nuclear population 
with 7 - 10^ (i.e., PeV energy /nucleon) necessary for 
matching the HEGRA data. For simplicity, for this pur- 
pose, we approximate the sum in Eq. (jl6p with a single 
species (take silicon). Setting the volume V = (4/3)7rr^ 
and the spectral index a = 2, we find from Eqs. ([T]) and 
(I17p the normalization constant 



N 



28 



1.0 X 10"^" ciR-^ TeV" 



(21) 
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E, [TeVJ 



FIG. 4: Energy-weighted 7-ray flux superposed over the 
HEGRA data. The curves are not actual fits to the spec- 
trum but are based on particular choices of the parameters 
which provide eyeball agreement with the data. The solid line 



corresponds to Si and the dashed one to ' 



N21 



2x 10" 



"cm ^TeV"^ and iVse 



10 



'Fe. We have used 
'i°cm-^TeV-^ In 



both cases, -E";^^ — 1.5 MeV and the spectral index a — 2. 

Note that the choice El^j^ — 1.5 MeV differs from the value 
1 MeV which was used as fiducial input for the discussion in 
the text. 



From Eq. (I14p this gives a required Si density at 1 PeV 



-{En = IPeV) = — -te^{En 



dEsi 



28 dE 



N 



IPeV) 



3.3 X 10"^^ cm-^ TcV"^ (22) 



of the Galaxy in order to provide a harder spectrum 
oc E~^. In order to ascertain whether this is possible, 
we compare the total energy in the nucleus population 
with the wind energy budget ~ 10'^^ erg/s. Integration of 
ENd.nsi{EN)/dEN, with the use of Eq. ((2T|) . gives for the 
volume of the HEGRA source a total energy ~ 10'*^ erg 
(modulo a logarithmic dependence on the source cutoff) . 
This energy is accumulated over the diffusion time of 
10'* yr = 3 X 10^^ seconds, so that the required power 
density is 2 orders of magnitude smaller than the kinetic 
energy budget of the entire association. 

In summary, we have presented in detail the criteria for 
the model to be viable. If future observations should lend 
credence to the model, significant insights into the accel- 
eration mechanism can be obtained, especially to provide 
clues with respect to the acceleration efficiency. In com- 
pleting the explanation of the HEGRA signal, there is 
of course one issue that needs to be addressed - the sig- 
nal was observed only in a 3 pc radius cell at the edge 
of the association. With the same angular radius, there 
are a total of 37 cells in the core of the association. The 
flux in each of the other cells is limited to be < 0.17 of 
TeV J2032-I-4130. At this point in our understanding we 
can provide only qualitative remarks. The obvious pos- 
sibility is an increased density of very hot OB stars in 
the cell of TeV J2032+4130 which provide more efficient 
trapping and accelerating conditions, as well as a hot- 
ter photon background. Indeed, a recent estimate [l^l 
indicates around 10 O stars in the region of the source, 
which is a factor of 3 larger than expected on the basis of 
a uniform population. (This strongly supports our previ- 
ous choice of R\-) More understanding of this can only 
come with more data and a full Monte Carlo simulation. 



It is of interest to compare this with the observed diffuse 
cosmic ray nuclear density dncR,si/dE in this energy re- 
gion. From fn\ . we find 

' IPeV) = — JcR.si 1 PeV 
dE c 

« 1.3 X 10"^^ cm"3 TeV~\ (23) 

a difference of 5 orders of magnitude. Combined with 
the observation that the diffuse spectrum has a power 
index of 3, this shows explicitly that (i) the required 
population density at PeV/nucleon must arise from ac- 
celeration of much more plentiful nuclei that are trapped 
at much lower energies (ii) the acceleration efficiency of 
the association must be considerably greater than that 



D. TeV Neutrino Emission 

As mentioned above, the interaction of nuclei with 
the background photons produces a beam of neutrons. 
The decay mean free path of a neutron is C7r„ = 
9.15 {En/10^ GeV) kpc, the lifetime being boosted from 
its rest-frame value r„ = 886 seconds to its lab value 
via 7 = En/mn- This means that for a source distance 
d ~ 1.7 kpc, practically all neutrons with E^ ~ 10^ GeV 
will decay en route to Earth, producing a flux of direc- 
tional antineutrinos. 

The basic formula that relates the neutron flux at the 
source {dFn / dEn) to the antineutrino flux observed at 
Earth {dF-/dET,) is Q 



dEn -j-priEn) 

atjn 



" dP , dcoa0v 



r 



S [Ett - £;„ ei7 (1 + cos 0T7)/mr. 



(24) 



The variables appearing in Eq. (I24p are the antineutrino and neutron energies in the lab (iJjj and £"„), the antineu- 
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trino angle with respect to the direction of the neutron 
momentum, in the neutron rest-frame (^ij), and the an- 
tineutrino energy in the neutron rest-frame (eu) . The last 
three variables are not observed by a laboratory neutrino- 
detector, and so are integrated over. The observable E- 
is held fixed. The delta-function relates the neutrino en- 
ergy in the lab to the three integration variables. The 
parameters appearing in Eq. (|24p are the neutron mass 
and rest-frame lifetime (m„ and t„). Finally, dP/dejj 
is the normalized probability that the decaying neutron 
produces a V with energy e^r in the neutron rest-frame. 
Note that the maximum v energy in the neutron rest 
frame is very nearly Q-p- = m„ — rup — nie — 0.78 MeV. 
Integration of Eq. ([24|) can be easily accomplished, es- 
pecially when two good approximations are applied. In 
the lab, the ratio of the maximum V energy to the neu- 
tron energy is Q— /to„ 10^^, and so the boosted v 
has a spectrum with E-jj e (0, 10^'^£'„) and an average 
energy eo = 0.48 MeV. The first approximation is to sim- 
plify the antineutrino spectrum from /3-decay and take a 
monochromatic spectrum with energy equal to the av- 
erage energy, eq. In such a way dP/de-jj = 5{e-cF — eo) 
and Eq. simplifies. The second approximation is to 
replace the neutron decay probability 1 — e"''™"/^"'^" 
with a step function Q{E^'^^ — En) at some energy 
ijmax _ 0{m„d/Tn) - (d/9.15 kpc) X lO'^ GeV. Com- 
bining these two approximations, one obtains 72| 



dF- fJln 

dEkr 2 eo 



dEn dFn 



(En) . (25) 



The neutron emission can be related to the HEGRA 
7-ray flux only if photo-de-excitation following photo- 
dissociation is unsuppressed. Thus, a lower bound for 
neutron emissivity (and resulting neutrino flux) may be 
obtained by finding this relation for the cases (Si/Fe) 
discussed in the previous section. An extremely impor- 
tant example where there is 100% suppression is the case 
of ^He, because neutron emission results in transitions to 
stable A = 3 daughter states [t^ . This possibility will be 
discussed more at the end of this section. Meanwhile, we 
turn to estimate the neutrino event rate associated with 
Si dissociation which can be expected at the IceCube de- 
tector, now under construction at the South Pole [t^ . 

For a transition in which an average of riis 7-rays are 
emitted during de-excitation, an average of ~ 1/2 neu- 
tron is emitted during the stripping (there is almost equal 
probability of emission of n and p). The conservation of 
Lorentz factor allows the relation 



7B' 



2« dF^ 



dE^ 



dE~i = 2 



dFn 
dEn 



dEn 



(26) 



Taking the derivative with respect to 7 gives the relation 



El 



dF^ 



dE^ 



2 rUn n28 



dFn 



dEn 



(27) 



Eji—jmj^ 



which leads to desired relation between photon and neu- 
tron fluxes 



dFn 

dEn 



(En) 



dF., 



2m„n28 dE^ 



(E^ = EnE'/nin) . (28) 



Substituting into Eq. (|25p . one finds the antineutrino flux 
associated to a given flux of photons 



dF^ 1 



2eo 

.^728, 



1 



dF 



2n^ dE^^^'' ~ ' 



(29) 

where a is the spectral index of the photon population. 
When referring to the 3 pc HEGRA cell, and taking 
2eo ~ ^7281 finally obtain the predicted associated 
antineutrino flux 



dF- _ 1 TYhegra / g- \ 



(30) 



which is valid in the energy window for which nuclei dis- 
integration takes place, i.e., F^ ^ (10^ — 10^) eo. 

At IceCube, the events are grouped as either muon 
tracks or showers. Tracks include muons resulting from 
both cosmic muons and from Charged Current (CC) in- 
teraction of muon neutrinos. The angular resolution 
for muon tracks « 0.7° ^1] allows a search window of 
1° X 1°. This corresponds to a search bin solid angle of 
Ar^ioxio « 3 X 10"'* sr. Since IceCube does not resolve 
the 3 pc HEGRA cell, the V flux will have contribution 
from all 37 cells in the core of the association. To es- 
timated this, we assume that the flux from each of the 
36 non-HEGRA cells is equal to the upper limit found in 
the direction of Cygnus X-3, namely 1/6 of the flux from 
TeV J2032-f 4130. The antineutrino flux from the Cygnus 
region can then be approximated (perhaps generously) as 



dF-.^ ^ 7 A^HEGRA / Ft, 



^ / \ ' J »tlKGK 

dEp 4 7128 



VTeV 



(31) 



To estimate the expected number of induced tracks 
from Cygnus 0B2 we adopt the semianalytical calcula- 
tion presented in Ref. [76| . 



A/''"' 
■'^Cyg 



(32) 
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where dFj^^/dE- is the flux, dace I is the diff'er- 
ential CC interaction cross section producing a muon of 
energy i?", is the number density of nucleons in the 
matter surrounding the detector, and t is the exposure 
time of the detector. After being produced, the muon 
propagates through the rock and ice surrounding the de- 
tector and loses energy. We denote by F{E^, i?^", I) the 
function that represents the probabihty of a muon pro- 
duced with energy E'^, arriving at the detector with en- 
ergy i?^", after travehng a distance I. The details of 
the detector are encoded in the effective area ^eff- We 
use the parametrization of the A^jj described in Ref. [t^ 
to simulate the response of the IceCube detector after 
events that are not neutrinos have been rejected (this is 
achieved by quality cuts referred to as "level 2" cuts in 
Ref. [ij)- The minimum track length cut is Imin = 300 m 
and we account for events with > 500 GcV. 

Although the flux of antineutrinos produced by Cygnus 
0B2 is pure Fe, because of neutrino oscillations, the an- 
tineutrinos observed at Earth will be distributed over all 
flavors, 



dFj, 



dEj, 



5 + A'-HS' ™ 



where a = e, fi, t denotes the neutrino flavor. The 95% 
confidence ranges for the probability differences. 



were derived elsewhere [77[ using the results of the up- 
to-date "iv oscillation analysis of solar, atmospheric, LBL 
and reactor data [t^. Substituting Eqs. and (|34l) 
into Eq. ([5^ we obtain the F^-induced tracks in IceCube 
from Cygnus 0B2 in i = 15 years of observation, A/^yg « 
7.5/n^. 

Showers are generated by neutrino collisions — 
Ve or Vt CC interactions, and all Neutral Current (NC) 
interactions — inside of or nearby the detector, and by 
muon bremsstrahlung radiation near the detector. For 
showers, the angular resolution is significantly worse than 
for muon tracks. In our analysis, we consider a shower 
search bin solid angle, Ar2iooxio°- Normally, a reduction 
of the muon bremsstrahlung background is effected by 
placing a cut of 40 TeV on the minimum reconstructed 
energy [79]. For Cygnus 0B2, this strong energy cut is 
not needed since this muon background is filtered by the 
Earth. Thus we account for all events with shower energy 
Esh > E^'^ = 1 TeV. The directionality requirement, 
however, implies that the effective volume for detection 
of showers is reduced to the instrumented volume of the 
detector, Voft — 1 km'', because of the small size of the 
showers (less than 200 m in radius) in this energy range. 
Following Ref. [TT'I we estimate the expected number of 
showers from Cygnus 0B2 as 



-0.082 I 

h^^^^ = -0.128ir'' 



0.1061°:°^^ 

" f^, , (34) 



where 



IV rsh 



. rsh.CC 



.rsh,NC 
'■^Cyg 



(35) 



(36) 



and 



dEL / dEh 



.rsh,NC , 1, 

^Cyg =tn^V.. ^ -dE^^-^'lEE:. 



(37) 



a— e./i,r 



Here, da^QjdE'— is the differential NC interaction cross 
section producing a secondary antineutrino of energy, EL. 
In writing Eqs. ([36]) and ([37]) we are assuming that for 
contained events the shower energy corresponds with the 
interacting Fe or antineutrino energy (iSgh = ^V) in 
a CC interaction, while for NC the shower energy cor- 
responds to the energy in the hadronic shower i?sh — 
Eir — E!jj = y Etj where y is the usual inelasticity pa- 
rameter in deep inelastic scattering. In total we expect 
•^Cyg ~ 5/rT28 from Cygnus OB in 15 years of observa- 
tion. 

We now turn to the estimate of the background of 



atmospheric neutrinos. For the "conventional" atmo- 
spheric neutrino fluxes arising from pion and kaon de- 
cays, we adopt the 3-dimensional scheme estimates of 
Ref. [sOl, which we extrapolate to match at higher ener- 
gies the 1-dimensional calculations by Volkova [8l|. We 
also incorporate "prompt" neutrinos from charm decay as 
calculated in Ref. 82]. We obtain the number of expected 
track and shower events from atmospheric neutrinos as 
in Eqs (1221), dSSl), and ^ with dF^™/dE^{E^) being 
the Vf, and u^^ atmospheric neutrino fluxes integrated over 
a sohd angle of 1° x 1° (for tracks) and 10° x 10° (for 
showers) width around the direction of the Cygnus 0B2 
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Primary quanta 


7 production mech. 


v production mech. 


Number ratio 


Energy ratio 


comments 


A 


A~* A* rTJ-y + A 


A -+ n — > 77g 
[E^ ~ 7.4 to] 


7 : ~ riA ■ 1/2 


E^ '^-fA 
E^ ~ Ef, 

~2-8 


observation of v's at IceCube 
depends on ''He abundance 


P 


p ^ 7r° ^ 2 7 


p TT^ Siy 


7 : '--J 1 : 3 




i^'s ARE seen at IceCube 


e- plasma 


synchrotron 
inverse Compton 


none 









TABLE I: Comparison of 7-ray and neutrino emission from A, p, and e primaries. Note that per 7-ray, an order of magnitude 
fewer neutrinos are expected from nuclei photodisentigration than from hadronic interactions followed by pion decays. Note 
also that the neutrino energy from the nuclei photo-disintegration is typically about one order of magnitude smaller than the 
7-ray energy. When the primaries are electrons, only 7-rays are produced, but not neutrinos. 



source 9 — 131.2°. We get an expected background of 
MX^yy = 14 and A^Itm = 47 in 15 years. Of the 47 
showers, 16 correspond to Ve CC interactions while 31 
correspond to NC interactions. 

It is clear from the preceding that there is no significant 
antineutrino signal resulting from dissociations which are 
accompanied by unsuppressed photo-de-excitation. This 
is in complete contrast with the pp case where for an 
spectral index of 2 one ex pec ts on av erag e a flux of 
half of that of 7-rays [H, [H, [H [H, [83| (see Sec. E]). 
Then, when considering the 37 cells of the Cygnus 0B2 
association one expects a Vfj^ event rate of 4.2 yr~^ with 
a background of 0.9 yr^^. The track signal in the pp 
case is about an order of magnitude larger than in the 
A* case. Because of production and oscillations the pp 
mechanism yields about one per 7-ray. The yield in 
the A* mechanism (for nis = 2) is about 0.25 Ve per 
7-ray which, after oscillation, yields 0.05 Vf^ per 7-ray. 
This ratio decreases by an additional factor of ~ 2 if 
comparison is made with integrated spectrum above 1 
TeV. 

We turn now to comment on the role of helium. Except 
for protons, helium nuclei dominate the cosmic ray spec- 
trum, with a population about 100 times larger than the 
heavy species [H, HI]. As mentioned above, the strip- 
ping of a nucleon from ^He leaves the residual A — i 
nuclides in their ground states [t^, so that there is no 
photon emission. However, stripping to "^He with emis- 
sion of a neutron will provide a yield of neutrinos. As 
can be seen from Fig. [31 the stripping rate is down by an 
order of magnitude from that for the heavy elements so 
that, with the larger population, the expected antineu- 
trino flux would be about a factor of 10 larger than our 
prediction from Eq. (j3ip . This enhancement has pow- 
erful consequences for the source energetics: the energy 
required to accelerate the entire *He population is very 
close to the allowed energy budget for the diffusion life- 
time. Therefore, should IceCube obtain a statistically 
significant signal that cannot be ascribed to pp interac- 
tion (because of observation of the TeV 7-ray suppression 
predicted by the model), then there could be a hint of 
extraordinary efficiency in the trapping an accelerating 
mechanisms in extremely hot and intense stellar associa- 



tions. 



V. PHOTO-DISINTEGRATION vs tt DECAY 

In the previous section we have described the mecha- 
nism of 7-ray production from nuclei de-excitation af- 
ter disintegration in the background photon field and 
showed that the HEGRA data could be explained in 
these terms. In addition, this mechanism can give rise 
to a neutrino flux after the stripped neutrons decay in 
flight. However, as we mentioned in the Introduction, 
there are two channels other than photo-disintegration 
that might contribute to 7-ray and neutrino production 
at the Cygnus region. These are photo-hadronic (A-7) 
and pure hadronic {A-p) interactions. In both cases, 7- 
rays (neutrinos) are produced after tt^ (tt^ and tt^) de- 
cays. 

As noted above, photo-meson production has a very 
high energy threshold, being only relevant for very high 
energetic beams or in very hot photon enviroments. Even 
in these extreme cases, the fact that this reaction turns on 
at so high energies implies that the photons and neutrinos 
from decaying pions are produced at very high energies 
too, well above the TeV range. Hence, in the following 
subsections we comment on the 7-ray and neutrino emis- 
sivities due to nucleus-proton collisions and compare it 
with those calculated for the photo-disintegration of nu- 
clei. But first, in Table U we show approximate estimates 
for the energy, number ratio and energy ratio of neu- 
trinos and 7-rays due to photo-disintegration and pion 
decay at production (we have also included the leptonic 
mechanism of 7-ray production for completeness). Now, 
we will compare these two hadronic mechanisms. 



A. TT spectrum 

The interaction of high energy nuclei with the cold 
ambient interstellar medium (ISM) gives rise to 7-rays 
through the decay of the produced neutral mesons. The 
tt'^ emissivity resulting from an isotropic distribution of 
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accelerated nuclei dn{EN)/ dE^ is given by [9C 



^ dn , . 

d(7A 



dE^o 



{Et^o,En) dEjy 



(38) 



where tt-h is the ISM number density, E*^{E^o) is the 
minimum energy per nucleon required to produced a pion 
with energy E^^a, and daA^E^o^ Eiq)/dET^a is the differen- 
tial cross section for the production of a pion with energy 
E^^o in the lab frame due to the collision of a nuclei A of 
energy per nucleon Em with a hydrogen atom at rest. 

Hence, an accurate knowledge of the differential cross 
section for pion production is necessary to calculate the 
7-ray emissivity from this channel. There have been sev- 
eral approaches and parameterizations in the literature, 
most of them based on the use of the isobaric [9l[ and 
scaling fo^, 'os'l models of the reaction or their combi- 
nation [94i . .95 . 96]. The (5-function approximation was 



considered in Ref. [97| and the inclusion of diffractive 
interactions and scaling violations in Ref. |98j| (see the 
Appendix of Ref. [9^ for a a comparison of different ap- 
proaches). A new recent parametrization is presented 
in Ref. [100l | based on simulations of proton-proton in- 
teractions from the SIBYLL event generator [lOl]. The 
isobaric model has been shown to work reasonably well 
at low energies [E <Z GeV), whereas the scaling model 
is more suitable at higher energies. Hence, here we will 
follow the scaling model with a parameterization of the 
differential cross section which is an_Ejv-independent ap- 
proximation of that given in Ref. [lOOj ] 



d(TA 

dE^o 



iE^o,E 



N) 



E^ 



(39) 



where x = Et^o/En and = A'^/^ctq, with cto ~ 34.6 
mb, which takes into account the scaling of the cross 
section with the atomic number 102J, and 



J 



/^o(x) ~ 8.18x^/2 



1 - 



1.33(1-2x1/^) 

1-h 1.33x1/2 (l_2;l/2)y \l-x^/^ ^ 1 + 1.33x1/2 (l_2;l/2 



(40) 



r 



which takes into account the high pion multiplicities at 
high energies. 

By using this form for the differential cross section and 
a power-law cosmic-ray spectrum, the ir^ emissivity can 
be written as 



where E^]i'^{E^) ^ E^ + mlo/iAE^). Hence, the 7-ray 
emissivity is given by 



]^P{E^)^Z^o^(a)Q^i{E^) 



(45) 



where 



rSiE^o)^ZAMa)QfiE"^) 



}f(EM)=<7^cn^i^iEM) 



dn 
dE^ 



(41) 



(42) 



with Zj^o^[a) — 2/a. 

On the other hand, the 7-ray emissivity due to nuclei 
photo-disintegration, Eq. ()16p . in the range of energies 
where Ra is approximately constant, can be written as 



and the spectrum- weighted moment of the inclusive cross 
section or so-called Z-factor is given by 



Q^'^iE^) ^ ZA-,{a) QfiE^ 



(46) 



Zati<> (a) = / x" ^ /^o (x) dx 
Jo 



(43) 



where, assuming dominance of a single component (sili- 
con), 



where, as usual, a is the spectral index of the cosmic-ray 
spectrum. 



B. TeV 7-rays 



Since isotropy is implied in Eq. psp . the 7-ray emis- 
sivity is obtained from the tt^ emissivity as 



^'iE,) = 2 



THdE^o (44) 



QtiE. 



N 



ZA^{a) 



TUN 



a-1 



(47) 
(48) 



Hence, the ratio of the 7-ray emissivity due to these 
two mechanisms is given by. 
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-) 

28 / 



3/4 



10-5 yrs" 



0.005 \ / Zj^^a{a) Z^o^{a) 



(49) 



Thus, for riH = 0.1 cm-'^ as quoted in the model of 
Ref. [48| , 7-rays from nuclei dissociation dominate those 
from pion decays produced in hadronic interactions by 
as much as an order of magnitude. Note however that 
in the model of Ref. [iMi the volume-integrated luminos- 



ity / Q^P{E^)dV cx n. 



-1/2 



Thus, within this model. 



a higher density would decrease the hadronic 7-ray flux 
coming from this region. As suggested in Ref. [48|, we 
have taken nu = 0.1 cm"^ which, for the concentration of 
stars in the Cygnus 0B2 region, corresponds to a volume 
of interaction equal to that of the region of the HEGRA 
signal. 

Hence, 7-rays produced by the de-excitation of nuclei 
after photo-dissociation in the photon background might 
well be the primary contribution to the TeV signal from 
the Cygnus 0B2 region. 



C. TeV neutrinos 

The interaction of high energy nuclei with protons of 
the ISM gives rise to a neutrino (and antineutrino) emis- 
sivity at production which, in an analogous manner as in 
the previous case, can be expressed as 

Q^P{E,J = 2 [Z^o^^ia) + Z^,^{a)] Q^!^{E,J (50) 

QiP{E,J = 2Z,,,^{a)Q^P{E,J (51) 

where the factor 2 takes into account the sum of both 
neutrino and antineutrino emissivities, and Z^^o^^^a), 
Zfuj^ (a) and Z^^^ (a) are the Z-factors corresponding to 
from pion and muon decay and Ve from muon decay, 
respectively. They are given by ^lOQ] 



Z^rO^^ia) 



(1-r) 



a-l 



(52) 



4 [3-2r-r° {3~2r + a^ar), ^^^^ 



a2 (1 - r)2 (a -f- 2) (a + 3) 

24 [a{l-r) -r{l-r")] 
(l-r)2(a + l)(a + 2) (a + 3) 



(54) 



where r — {m^/m^Y' — 0.467. Hence, after taking into 
account oscillations, the ratio of muon and electron neu- 
trinos plus antineutrinos to 7-rays of the same energy 
for a £ (1.5,3) is approximately given by (and in good 
agreement with Refs. [13, H Ha, S 113 ) 



R^^, (a) 



1.05 - 0.28a 



1.08 - 0.28 a 



(55) 



(56) 



Thus, the predicted ratios for v^w^-.Vt-:^ at a common en- 
ergy are '^1:1:1:1; here we have included v^^-Vt equilibra- 
tion which gives the added information that i/^:i'^=l:l. 

Note that at the high energies we are interested in, the 
ratio of the number of neutrinos to that of 7-rays at the 
same energy, is smaller than the expected 3:1. The reason 
being that the average neutrino energy is about a quarter 
of that of the parent pion, while in the case of the photon 
it is half. Hence there are more very low energy neutrinos 
than very low energy photons, which compensates what 
happens at higher energies, so that the ratio of the total 
number of neutrinos to that of photons is indeed 3:1. 

On the other hand, for the nuclei disintegration case, 
from Eq. (|29p . the neutrino emissivity at production can 
be written as 

QtS^vi) = R%{a) Z^.{a) Qf (i?^J (57) 

where the Z-factor associated to the antineutrino emis- 
sion from neutron decay and the ratio of neutrons to 
7-rays of the same energy, R'^J^ , are given by 



Znu{.a) 



2n28 



/2eo 



. ^728 , 



a-l 



(58) 



(59) 



Hence, the ratios of neutrinos (after taken into account 
oscillations) due to the Ap and photo-dissociation mech- 
anisms are given by 
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IQ-Syrs-i 
Ra 

R^/%{2) ZU^) Za^{2) 
i?f-^(a)Z„,(a)Z^^(a) 




Zatto{2) Z^o^{2) 



(60) 



-RAp/dis(") 




(61) 



i.e., (and z^t) is suppressed by the factor 4.2 in photo- 
dissociation relative to Ap, and i/^ is suppressed by 1.8. 
We note that the neutrino flavor results here agree with 
the well-known neutrino flavor ratios ~ 5:2:2 and 1:1:1 
from astrophysical neutron and pion decay. 

Therefore, if the TeV 7-ray signal from the Cygnus 
0B2 region can be explained by de-excitation of daughter 
nuclei after photo-dissociation, the neutrino population 
will be, in general, dominated by pion decay. However, 
with a significant population of ^He, this may not be the 
case: through photo-dissociation "^He only contributes to 
the i^-flux without affecting the photon signal. Pho- 
ton production from '^He p interactions are negligible, 
see Eq. P^j) . 



VI. EeV COSMIC RAY EXCESS? 

A seemingly different, but in fact closely related sub- 
ject we will discuss in this section is the intriguing 
anisotropy that has been observed by several experiments 
in the energy range near 10^ GeV = EeV. Analyses of 
the data collected by the Akeno Giant Air Shower Array 
(AGASA) seem to indicate a 4(t excess of events from 
the direction of the Galactic pla ne in the narrow energy 
window 10*-9 GeV to lO^-^ GeV ^0^. A Galactic plane 
enhancement at the 3.2a level between lO* '^ GeV and 
10^'^ GeV has a lso been reported by the Fly's Eye Gol- 
laboration [lo3|. The anisotropy in the AGASA data 
sample is dominated by hot spots near the Galactic cen- 
ter and the Gygnus region. However, based on a data 
sample larger than any previous experiment, the Pierre 
Auger GoUaboration has reported no evidence for such 
an excess from the direction of the Galactic center jlOSl ] . 

The complete isotropy below about 10^'^ GeV revealed 
by KASGADE data 106] vitiate direction-preserving 
photons as primaries. Therefore, the excess from the 
Galactic plane is very suggestive of neutrons as candi- 
date primaries, because the directional signal requires 
relatively-stable neutral primaries, and time-dilated neu- 
trons can reach the Earth from typical Galactic distances 



when the neutron energ y exceeds 10^ GeV. The analysis 
of Fly's Eye data [104| implies that if neutrons arc the 
carriers of the anisotropy, there needs to be some contri- 
bution from at least one source closer than ~ 2 kpc. As 
we mentioned in the Introduction, the full Fly's Eye data 
include a directional signal from the Cygnus region which 
was somewhat lost in an unsuccessful attempt [T2j to re- 
late it to 7-ray emission from Cygnus X-3. This signal 
was also observed with the Akeno air shower array (see 
Fig. [1]) 13]. As shown in Fig. [2l Cygnus 0B2 overlaps 
along our line of sight with Gygnus X-3, which lies about 
8 kpc farther away than the stellar association. 

In the cont ext of Cygnus X-3, we note that the 1989 
analysis jlOTi ] of data collected in the Haverah Park ex- 
periment was consistent with the absence of anisotropy in 
the direction of that source, in the energy bin of interest. 
There was indeed an excess, with a Poissonian probabil- 
ity of 0.013, but at energies above 10^-^ GeV. How ever, a 
more recent analysis by the group at Leeds [l08| ]. based 
on improved shower simulations, more selective trigger 
criteria, and a reduction of the z enith acceptance from 
60° to 45°, reveals that previous [l09l ] energy estimates 
need to be reduced by about 30%. It is then suggestive 
that the excess may overlap the region of interest. One 
can encourage an event-by-event analysis to see whether 
this is indeed the case; in the meantime one may re- 
serve judgment with respect to presence or absence of 
anisotropy in the Haverah Park data. 

More recently, the HiRes Collaboration has carried out 
a systematic search for point sources in the Northern 
sky. In an initial series of papers, they report nosig- 
nificant excess for energies greater than 10^'^ GeV [llOl |. 
10^" GeV [nil], and 10^°-^ GeV Tiaj. All these energy 
thresholds exceed the range above 10^'^ GeV where the 
Akeno and Fly's Eye anisotropics were reported. Mo- 
tivated by the present study, we have asked the HiRes 
Collaboration to extend their analysis down to the en- 
ergy threshold of interest, in the direction of Cygnus. 
The result is that there is no significant excess above 
background for the same angular and energy cuts used 
by Akeno and Fly's Eye, this with a larger data sample 
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FIG. 5: Photo-disintegration rate of ^"^Fe on the T = 50 K 
Cygnus OB2 blackbody radiation as given by Eq. ([8}. 



than available to either of these previous studies [lie 
Additional scans were performed for nearby energy and 
angular cuts, with the same result of no significant excess. 
If the previously-reported anisotropics were not statisti- 
cal fluctuations, then they can be associated with bursts 
over periods of less than 15 years. This possibility can 
be monitored with data to be collect ed at the Northern 
site of the Pierre Auger Observatory jll4| . 

In what follows, we will proceed on the ass ump tion 
that the anisotropy reported by the Fly's Eye [I4I and 
Akeno [l^ collaborations originated in a burst taking 
place in Cygnus 0B2 during a 10-yr period, and assess 
the implications for source energetics. 

For nuclei with energies EeV/nucleon, the dissocia- 
tion proceeds via collisions with far infrared thermal 
photons populating molecular clouds with temperature 
~ 50 K 115]. The photo-disintegration rate is shown in 
Fig. [5l giving an interaction time of about 100 yr, which 
allows a perturbative treatment during a 10 yr burst. We 
can then use Eq. ([9])with the left hand side set equal to 
the neutron luminosity for energies above Ei — 0.5 EeV. 
Specifically, 



dN 



dF„. 



El 



dEn 



dE„ 



7 X 10^'' neutrons s ^ , 



where 



dFn 



I El 



dE„ 



dEn ~ 2 X 10 



cm-2 



(62) 



(63) 



is the flux observed at Earth [l3, llJli and the factor 
gd/4.6 kpc corrects for neutron /3-decay. Now, taking 



Cyg 
'56 



10 ^ yr -'^ we obtain 



dn 



Fe 



dE 



dEN = 5.6 X 10"^° cmT^ . 



N 



(64) 



Assuming a nucleus power law spectrum dEpo/dE = 
iVse E~'^, Eq. ([64|) will determine the normalization con- 
stant. Then, the energy density can be obtained via 



El 



E^dE, 
dE 



(65) 



and so straightforward calculation shows that the total 
energy AttR^E/S = 5 x 10'^^ erg is well within the kinetic 
energy budget of a 10 yr burst. 



VII. CONCLUSIONS 

Over the last few decades, the Cygnus region has been 
the focus of extensive study with a view to understand- 
ing astrophysical phenomena in very energetic environ- 
ments. Interesting observations have been made, but 
have not held up at statistically significant levels. More 
recently, improved efficiency and higher resolution have 
provided signals at discovery level. Perhaps, the most 
precise of such signals is the (1999-2002) HEGRA detec- 
tion of the unidentified TeV J2032-I-4130 source, with an 
average flux ^ 3% of the Crab Nebula [lB|. This bright 
"hot spot" , at the edge of the very active Cygnus 0B2 
star association, has been recently conflrmed through a 
re-analysis of "old" data collected with the Whipple Ob- 
servatory [1^. However, the average flux (~ 12% of the 
Crab) detected during 1989-1990 is well above that re- 
ported by the HEGRA Collaboration. Neither the Whip- 
ple nor the HEGRA experiments see any evidence for 
variability within their individual databases. The large 
differences between the flux levels cannot be explained as 
errors in estimation of the sensitivity of the experiments 
because they have been calibrated by the simultaneous 
observations of other TeV sources. Certainly more data 
are needed to resolve this issue. Happily, data are cur- 
rently flowing from the Milagro telescope [ll6| | and the 
Tibet Air Shower Array ^llTj experiments at a significant 
rapid rate. 

In this paper, we have focussed on the unidentified TeV 
7-ray source observed by the HEGRA Collaboration. Be- 
cause of the absence of X-ray or radio counterparts, it 
has been difficult to provide a compelling explanation 
for the origin of the 7-rays from the HEGRA source. 
This has motivated us to present an alternative model, 
which has some well-defined predictions for future obser- 
vations, both in 7-ray facilities and at neutrino telescopes 
such as IceCube. The model proposes that the observed 
TeV 7-rays are the result of Lorentz-boosted MeV 7-rays 
emitted on the de-excitation of daughter nuclei follow- 
ing collisions of PeV nuclei with a hot ultraviolet photon 
background. Our results include the following: 
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(1) There is a specific prediction of a suppression on 
the 7-ray spectrum in the region < 1 TeV - this because 
of the need to achieve the Giant Dipole Resonance for 
disintegration through coUision with ~ few eV photons. 

(2) The flux of 7-rays produced through nuclei de- 
excitation process dominates by an order of magnitude 
that resulting from tt*^ production and decay in collisions 
of nuclei with the ambient gas background. 

(3) The presence of a significant neutrino signal from 
decay of neutrons produced in the course of the photo- 
disintegration is strongly contingent on the abundance 
of ^He in the source, mirroring its observed prevalence in 
cosmic rays. This compensates by an order of magnitude 
over its reduced photo-disintegration rate, and results in 
a Vfj^ signal at IceCube (after oscillation from Ve) com- 
parable to the v^j, signal from hadronic interactions (as 
noted previously, this is sufhcient to establish a 15(t dis- 
covery in 15 years of observation.) The ability of IccCube 
to measure flavo r ratios can distinguish between these 
possibilities [llSj . 

(4) For this model to succeed, a high efhciency is re- 
quired for accelerating a low energy population of nuclei 
at the source to PeV energies. Thus, evidence for the 
model (such as a substantial suppression in the 7-ray 
spectrum below 1 TeV) could provide some insight 
into source dynamics. 

(5) It is noteworthy that the contribution to the dif- 
fuse 7-ray flux component resulting from the interac- 
tion of extreme-energy cosmic ray nuclei with the back- 
ground radiation fields permeating the univ erse [119] is 
well below the observed EGRET data [113 (for details 
see Ref. 

In summary, the availability of new data in coming 
years in 7-ray, X-ray, neutrino and cosmic-ray sector ex- 
periments affords an outstanding opportunity for study 
of high energy processes in extreme environments. Care- 



ful exploration of model implications will provide an im- 
portant complementary role to these observations. 
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Note added 

In the accompanying Letter [T3| we have adopted a more 
descriptive notation for some variables. Throughout this 
paper we use the traditional notation, which shortens 
the text length of the equations. The equivalences are: 
7 = r^, e' = r^e, e'^ = cGDR, = gdxn^ ^ 

e^AB^ E = ^LAB^ r - PGDR, and (Tq = 
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